Coincident count rates in absorbing dielectric media 
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A study of the effects of absorption on the nonUnear process of parametric down conversion is 
presented. Absorption within the nonhnear medium is accounted for by employing the framework 
of macroscopic QED and the Green tensor quantization of the electromagnetic field. An effective 
interaction Hamiltonian, which describes the nonlinear interaction of the electric field and the linear 
noise polarization field, is used to derive the quantum state of the light leaving a nonlinear crystal. 
The signal and idler modes of this quantum state are found to be a superpositions of the electric 
and noise polarization fields. Using this state, the expression for the coincident count rates for both 
Type I and Type II conversion are found. The nonlinear interaction with the noise polarization 
field were shown to cause an increase in the rate on the order of 10 for absorption of 10% per 
cm. This astonishingly small effect is found to be negligible compared to the decay caused by linear 
absorption of the propagating modes. From the expressions for the biphoton amplitude it can be 
seen the maximally entangled states can still be produced even in the presence of strong absorption. 
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I. INTRODUCTION 

The nonlinear phenomenon of parametric down conver- 
sion (PDC) is ubiquitous across all areas of quantum op- 
tics. The strongly correlated photon pairs generated by 
this process have wide applicability in a disparate range 
of areas; from fundamental tests of quantum mechanics 
[iHQ to quantum cryptography and quantum informa- 
tion processing [BHS]. As a result, an extensive body of 
literature exists on this process (for a necessarily incom- 
plete selection, see Refs. [9HT8]). However, despite this 
continued interest, many features of PDC have not been 
studied in detail. As the focus is usually on the interact- 
ing electromagnetic fields, it is common to neglect some 
of the more complicated aspect of the nonlinear medium. 
As a result, absorption is usually ignored. 

Absorption, however, plays a key role in all physical 
(response) processes as it is an essential requirement of 
causality. Methods for consistently quantizing the elec- 
tromagnetic field in absorbing linear electric and mag- 
netic materials have been developed (for a recent review, 
see Ref. [E]), and recently this formalism has been ex- 
tended to include nonlinear processes. In Ref. [20 the 
authors derive an effective interaction Hamiltonian to 
characterize PDC in absorbing media. The result showed 
the appearance of extra noise field interactions, with the 
pump photon now having the ability to convert not only 
to electric field modes, but also to noise polarization 
modes or a combination of the two. Furthermore, the 
Hamiltonian leads to the consistent inclusion of linear 
absorption as well. Hence, a full description of the effect 
of absorption can be obtained using this formalism. 
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In this article we will consider the effect of these ex- 
tra features on the coincident count rate for Type I and 
Type II down conversion processes. In Sec. [n| by study- 
ing the evolution of the quantum state of the signal and 
idler modes subject to an effective nonlinear interaction 
Hamiltonian, we derive the state vector for the outgoing 
modes. In Sec. [1111 the state vector will be used to calcu- 
late the coincident count rates for both types of crystal 
configuration. Discussion and some concluding remarks 
can be found in Sees. [TVl and [Vl 



II. THE EFFECTIVE INTERACTION 
HAMILTONIAN AND THE STATE VECTOR 

We begin by briefly reviewing the concepts of field 
quantization in absorbing (linear) dielectric media. The 
electromagnetic field is expanded in terms of the classical 
Green tensor for the Helmholtz operator [19 . The macro- 
scopic fields can then be expressed in terms of bosonic 
operators that describe collective excitations of the elec- 
tromagnetic field and the absorbing matter. As a result 
the electric field operator becomes 

CO 

E(r) = J dujF^{r,uj)^ h.c. (1) 



with frequency components 

E(r,cj) = ^ / (i^5G(r,s,u;) • Pn(s,u;) 
= n/^^ / ^/?^MG{r, s, uj) • f (s, uj), (2) 

V SoTT J 

and e'\s^uj) being the imaginary part of the complex 
permittivity £{s^uj) = £'{s^uj) + ie"{s^uj). The Green 
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tensor (or dyadic Green function) G(r, s, cj) is the unique 
solution of the Helmholtz equation for a point source 

V X V X G(r,s,cc;) ^£:(r, cj)G(r, s, cj) = ^(r-s). (3) 

The frequency components of the noise polarization field 



P^{r,oJ)=i^J^e"{r,oJ)f{r,co) (4) 

drive the electromagnetic field and are due to absorp- 
tion processes inside the dielectric medium. The bosonic 
operators f(s,a;) and ^(s,^^) obey the commutation re- 
lation 

f(r,a;),ft(s,a;Oj = S{r - s)S{uj - ou'). (5) 
Finally, the bilinear Hamiltonian 



generates the time-dependent Maxwell equations from 
Heisenberg's equations of motion for the electromagnetic 
field operators. 

Let us now consider the setup schematically depicted 
in Fig. [I] A x^^^ crystal of length L (region 2), sur- 
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FIG. 1: A x^'^^ crystal of length L is pumped from the left 
by a continuous, single- frequency classical plane wave. The 
resulting signal and idler photons, with wave vectors qs and 
Qi, propagate out of the opposing face of the crystal and are 
detected in coincidence by detectors at rds and Tdi- 

rounded by free space (regions 1 and 3), is pumped from 
the left by a continuous, single- frequency classical plane 
wave with frequency Up. The classical pump field ap- 
proximation asserts that the pump field is of sufficiently 
high intensity such that it is negligibly depleted by the 
nonlinear interaction. Outside the crystal the pump has 
wave vector magnitude qp = Up/c. Inside the crystal the 
pump has wave vector magnitude kp = n(YA^ojp)ujp/c^ 
where n{vA^oj) is the refractive index of the medium. 
Owing to absorption, the refractive index, and hence the 
wave vector inside the crystal, is described by a com- 
plex number. The coordinate system is taken such that 
the pump beam propagates along the z-axis. The pump 
beam enters the crystal at the planar boundary between 



regions 1 and 2 located at z = —L/2. The interface 
plane wave (amplitude) refiection and transmission co- 
efficients are r^^M ^^em^ respectively. The pump 
beam interacts with the medium to produce signal and 
idler modes with frequencies uOg^i- The magnitudes of the 
wave vectors of these modes inside and outside the crystal 
are ks^i = n(YA,^s,i)^s,i/ c and qs^i = ujs^i/c, respectively. 
The signal and idler modes propagate out of the crystal at 
the planar boundary between regions 2 and 3 located at 
z = L/2. The (amplitude) reffection and transmission co- 
efficients for this interface are t'^^/tm ^'te/tm (^o^^ 
that the outgoing modes are no longer plane waves). For 
simplicity the crystal is assumed to be infinitely extended 
in the x — y plane. The outgoing modes are detected by 
coincidence counters located at Vds and Vdi. 

In Ref. ^20j we have shown that the effective Hamil- 
tonian that describes PDC in absorbing media can be 
written as 



Hf = J d^r J (icjfza;f'^(r,a;) •f(r,a;), (6) ^int = 



^^^^ Xaj^(^s,^i;^p) 



^i,(3 A , , ^) + >C [e* (r A , CJi )] PjJ -^^ {ta ,^i,t) 

Ep^^{rA,ujp,t) + jC[£{rA,ujp)]pNp,-f{rA,^p,t) 



h.c. 

(7) 



The greek vector indices, which refer to the polarization 
of the mode, run over the three Cartesian coordinates. 
Over these indices summation convention is implied. The 
function C[e{rA^oj)]^ given by 



C[£{rA,uj) 



2 e{rA,oj) - 1 
9^0 ^(i-A,^) 



(8) 



is the local-ffeld correction factor for the noise polariza- 
tion ffeld, and e{rA^oj) = n{rA,oj)'^ is the complex linear 
permittivity within the crystal. The spatial integral is 
taken over the volume of the crystal. 

The quantum state of light at the output face of the 
crystal can be found by evaluating the time evolution of 
the state vector subject to the interaction Hamiltonian in 
Eq. Formally, to ffrst order in perturbation theory, 
the long-time limit, steady- state quantum state vector is 
given by 

^ = 1^(0)) + ^ pJtH{t)\m) + o{Hi{tf). (9) 

The ffrst order truncation of the series implies that prob- 
ability of interaction is sufficiently low as to make the 
likelihood of a two consecutive down conversion processes 
negligibly small. 

The classical plane-wave pump ffeld can be found from 
the Green tensor description of a sheet current source at 
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infinity and is given by 

X le 



(10) 



written in a form familiar from classical wave optics. 
Here, ^ is the pump polarization vector and Zp is the 
distance from the pump source to the input face of the 
crystal. The Fresnel coefficients for reflection and trans- 
mission of plane waves at the boundaries are given by 



21/23 
^TEM 



(^) 



n{uo) — 1 



t 



12 

TEM 



(11) 



nH + l' nM + 1 

Multiple scattering within the crystal is accounted for by 

r(c^) = [1 - (12) 



TEM[ 



All these coefficients can be found in the far-field expan- 
sion of the generalized Fresnel and multiple scattering 
coefficients of the dyadic Green tensor for layered media 
(see Appendix [A|). 



The first term in the square brackets of Eq. ( 10 ) corre- 



sponds to photons that propagate directly to the interac- 
tion point. The second term corresponds to photons that 
are refiected at the back face of the crystal (see Fig. [2| . 
The crystal medium is taken to be in its ground state 
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FIG. 2: The first term in the square brackets of Eq. ( 10 ) cor- 



responds to propagation path (a), the second term to propa- 
gation path (b). Scattering off two or more interfaces is taken 
into account by the multiple scattering coefficient Mtem{(jo). 

and hence the (expectation value of the) noise polariza- 
tion field at the pump frequency vanishes. Furthermore, 
it is assumed that both the electric and noise polariza- 
tion fields at the signal and idler frequencies consist of a 
slowly varying amplitude and a fast oscillation 



^(s,i),a(rA,^s,i,^) = ^(s,i),a(rA,^s,i)e 
^N(s,i),a(l*A,^s,i,^) = ^N(s,i),a(l*A,^s,i)e' 



-iujs,it 



(13) 
(14) 



After combining Eqs. ([10|, ([13|), and ([14|) with the ef- 
fective Hamiltonian_|7|) , the resulting expression is sub- 
stituted into Eq. rfoFT^ The initial state for both the 



signal and idler modes are taken to be the vacuum 
(|?/^(0)) = |0)s|0)i). Finally, performing the integral over 
t results in a J-function. This enforces the energy con- 
servation condition ujp = ujs -\- oji. One should note that, 
as the frequency of a propagating mode is unchanged by 
a change of medium, this condition is true both inside 
and outside the crystal. Thus, one finds the steady-state 
quantum state vector at the crystal output as 



\^) =Ar|0),|0),- 



ih 



TEM 



{uJp)Mtem{(^p) 



[ dV^^,^(a;„a;,)[e^'^'" +e-^^-(^--^)4^^M(c^^)] 
Jv 



(15) 



where A/" is a normalization factor that is very close 
to one. Note also that a contracted notation for the 
nonlinear susceptibility has been used: 2daf3{cjOs,uJi) = 

Xalji^s.OJi'.UJpYep^^. 



III. COINCIDENT COUNT RATE 

Given a quantum state the count rate at the de- 
tectors located at r^^^ and Vdi is given by [TOl [TTj 

R= {ilj\El^^{rds,(^ds)El^{rdi,uJdi) 

X EsA^ds,^ds)Ei^i^{rdi,uJdi)\i^)' (16) 

Owing to the number of field operators in |?/^), one can 
rewrite the Eq. (16) as 



R = AUA^n 



(17) 



with 



^A/i = i{0\s{0\Es^x{^ds,^ds)Ei^,_,{rdi,uJdi)W (18) 

commonly being referred to as the biphoton amplitude. 

Using Eq. (15) to expand the state vector in Eq. (18) 
gives 



ih 



•,{0\Es,x{'rds,i^ds) 

El^{rA,uj,) + C[e*{rA,iOs)]PLj^A,^s)] \0)s 
Elp{rA,oJi) + C[e*{rA,uji)]P^,^l^{rA,uji)\ \0)i. (19) 
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The matrix elements can be evaluated using the commu- 
tation relations for the electromagnetic field operators 
given in Appendix [BJ One finds that the only contribut- 
ing terms are those for which uJds = and ujdi = oji. 



Upon applying the identity ImG = (G-G*)/2i, Eq. ([19) 
becomes 



,4 t^WK)MT£;MK)e'^^^^e^^2 

kpZA I ^— ik-p(zA~ L) 23 



with 



(20) 



(21) 



Let us consider the z components of the wave vec- 
tors of the three modes [c.f. Eq. (A3) in Appendix [A] . 
One can see that the four terms in Eq. (20) refer to 



four different types of directional scattering of the sig- 
nal and idler modes. The term proportional to GxaGfj^^ 
describes forward scattering where both of the outgoing 
modes leave in the same direction as the pump beam. 
The cross terms and conjugate term correspond to in- 
teractions where back-scattering of one or both of the 
outgoing modes occurs (see Fig. [3|. 

Closer inspection reveals that this wave vector depen- 
dence is contained in a number of exponential terms of 
the form e^'^fi.^^)^ ^ In general, for long bulk crystals, the 
dominant terms will be those for which f{kz) ^ 0, when 
the crystal is said to be longitudinally phase matched. 
Away from this condition the contribution to the total 
amplitude is proportional to s\ii(z[f{kz)\. Hence, even for 
relatively small deviations from perfect phase matching, 
the contribution to the total amplitude falls sharply. In 
order to conserve energy and be close to perfect phase 
matching, h 

thermore, since \kp\> \kz^s\-, \ kz 
comes from the terms where f{kz) ^ kp — kz^i — kz^s- All 
other terms lie far from phase matching and hence are 
negligible. One can see that the terms that are closest to 
this phase matching condition are contained in the term 
proportional to GxaG^^ in Eq. (20). Thus, the cross and 
conjugate terms can be neglected. In that way, one ex- 
tracts the relevant contributions to the nonlinear process 
in terms of Green functions, as a generalization but in the 
same spirit as the free-space analysis using plane waves. 

It is worth noting that, although negligible for long 
bulk crystals, the cross and conjugate terms may give 
comparable contributions in layered media [211 ES] . How- 
ever, the Green tensor method is well adjusted for the 
study of such systems as the Green tensor for layered 
media is well known and hence these effects can easily be 
investigated. 



kz^i and kz^s must have the same sign. Fur- 
the main contribution 
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FIG. 3: The interaction shown in (a) is the forward-scattering 
term given by the dominant term oc G\o,G^p. Interactions 
(b),(c) and (d) are described by the terms G\c,G^p^ GxaG*^^ 
and GJc^G*^, respectively. These involve back-scattering of 
one or more of the outgoing modes and can be neglected. 



Further to this, we shall also assume that the lin- 
ear permittivity is constant across the crystal, hence 
Aa{r,uj) = Aa{u;) in the range of integration. Thus, 
one can take these prefactors out of the spatial integral. 
As a result, Eq. ([2Q|) becomes 



A 



TTSo 



X / o^VAC^«;3(c^.,c^z)[e^'^"^+e"^'^^"^"'^V|.^^MK)] 
Jv 

X A*(cJs)A*(cJi)GAa(rds,rA,cJs)G^/3(rdi,rA,cJi) (22) 

which serves as the starting point for the special cases 
considered below. 

A. Type I down conversion 



The Green tensors appearing in Eq. (22) can be ex- 



panded in terms of vector wave functions as given in Ap- 
pendix |Aj As the crystal is infinitely extended in the 
(x, ?/)-plane, the integral over the transverse components 
of is proportional to a (^-function. Integrating over it 
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leaves one transverse -integral and the spatial integral 
in the z-direction over the length of the crystal. 

If the crystal is cut for Type I conversion the polar- 
ization unit vectors of the outgoing modes are parallel. 
Hence the nonlinear susceptibility has the form 



XX 
yy 



(23) 



ai3 



Here, x and y are the unit vectors for the polarization 
state of the signal and idler photons and are associated 
with the tensor indices a and p. The contraction of the 
TE and TM vector wave functions with such a suscep- 
tibility are given in Appendix [Cj In this geometry the 
cross-terms between A^(k^) and N(k±) cancel. 

Performing the z-integral over the length of the crystal 
and retaining only those terms that are close to perfect 
phase matching one finds 



(I) 



(24) 



with 



(Pk^ 1 



X e 



(27r)2 k\ k, 

^±-irds,±-rdi,. 



sin A^^ 



Nx (kx ) {-k±)XTM,TM {^s , ) 



ka ki 



(25) 



Here, 



X MTEM{^p)M^{uJs)M^^{uJi) 

X (1 + rfEM{^p)rl\oJs)rl}{uj,)e' 



(26) 



with cr, a' TE^ TM as required. Zds and Zdi are the 
perpendicular distances from the output face of the crys- 
tal to the detectors. The longitudinal phase mismatch 
and longitudinal phase sum are 



k, 

kr. 



(27) 
(28) 



with the z-component and magnitude of the wave vector, 
respectively, being kz = k = n{uj)uj/c inside 

the crystal and Qz = — k\^ q = uP' jc? \xi the vacuum. 



For a complete solution, the integral I^J} in Eq. (|25) 



A/i 



would have to be solved numerically. However, in cer- 
tain limits it is possible to obtain an analytical solution. 
For example, let us consider the degenerate (cj^ = cj^ = 
bj^l2 = uj, ks = ki = k), collinear {rds,^ = ^di,±) case. 



Expanding the tensor product of the TE and TM vector 
wave functions as shown in Appendix [C) and performing 
the integral over the angular component of k^ , one finds 
that terms with odd powers of kx or ky vanish. Hence, 



'A/i 



dk± k± 
47r /c? 



sin A A: f 
Akk 



Ia/x 



X 



TE,TE 



(^s,^i) 



X 



TM,TM 



(^s,^i) 



(29) 



where Ia/x are the components of the 2-by-2 unit matrix. 

An asymptotic approximation to the /cx -integral can 
be found in the far-field limit where one assumes that 
the distances from the crystal face to the detectors are 
large compared to the wavelength of the detected modes, 
zuj/c 1. As, on mere practical grounds, this is al- 
ways the case, this limit provides a good approximation 
for the count rate. However, one should note that for 
large crystal-detector distances the propagation of the 
outgoing modes is effectively parallel. Hence, in the far- 
field limit the count rate is asymptotic to the collinear 
case. For non-collinear propagation one cannot use this 
approximation and alternative methods for solving the 
integral must be found. 

In the far-field limit one can use the method of station- 
ary phase to evaluate the /c^-integral. At the stationary 
phase point at /cx = one finds that kz = k = n{uj)uj/c. 
Thus, the Fresnel coefficients become 

21/23/ N 21/23/ N ^(^) - 1 /on\ 

^TE (^) = -^TM (^) = ^^jy^' (30) 

.23 / X _ ,23 / X _ M^) 



l-r2i(w)r23(w)e2»(-)^^]"\ 

(32) 



and hence 



Xte,te{^,^) = Xtm,tm{^,^) = X^{uj,uj), (33) 
Xte,tm{^,^) = Xtm,te{^,^) = X-{uj,uj). (34) 



Furthermore, 



Ak = ^[n{2uj) - n(a;)], 
^k = ^[n{2uj) ^ n{uj)]. 



(35) 
(36) 



The leading-order contribution to the integral gives 



p'i^izds+Zdi) 



27r cjn(cj)2 {zds + Zdi) 
sin A A: f 
Akk 



e^^^^X+(c^,c^)lAM (37) 
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Substituting Eq. back into Eq. (24), one finds that with 



the biphoton amphtude for Type I conversion is given by 



sin Afcf 



(Zds + Zdi) 

e^^'=*X+(w,w)lA^, (38) 



with the corresponding coincident count rate 



\Epd{L0,L0)\^^ 



X \X+{oj,oj)\^ 



sinAfcf 



• (39) 



MI) ^ / « K± 



sin Afc| 



2 A;^ fcj, Ma (ki ) ( - k_L ) Xtb ,T B K , ) 



+ i{kl-kl) (^^^ M,(k^)iV^(- 



/Co /Co 



(42) 



with X^^/, AA: and SA: defined in Eq. ([26|), Eq. ([27|) and 
Eq. (I28|). 



For the degenerate, cohinear case, after expanding the 
tensor product of the TE and TM vector wave functions 
(Appendix [C|), and performing the angular integral, one 
finds 



B. Type II down conversion 



Similarly, for Type II down conversion the Green ten- 
sors in Eq. (22) are expanded in terms of vector wave 
functions as shown in Appendix |A] and the integral over 
the transverse components of r^^ is performed. In Type II 
conversion the polarization unit vectors of the outgoing 
modes are perpendicular. Hence the nonlinear suscepti- 
bility has the form 



xy 
yx 



(40) 



a(3 



As before, x and y are the polarization unit vectors asso- 
ciated with the tensor indices a and f3. The contraction 
of the TE and TM vector wave functions are again given 
in Appendix [C] One should note that for this geometry 
the cross-terms between the M(k±) and N(k±) TE and 
TM vector wave function do contribute. 

The z-integral is performed as before, keeping only 
those terms that are close to perfect phase matching. 
Thus, the biphoton amplitude becomes 



HE, 



f ,2, ,2 



(11) 



d{uJs,uj,)A%ujs)A%uj,)ll^'^ 



(41) 



All) _ f dk±k. 



87rfc? 



sinAfcf 



^i^k^ ^iqz {zds+Zdi) J 



X/j, 



(43) 



with the leading-order contribution to the far-field limit 
given by 



ai^{Zds+Zdi) 



sin A A: f 



47rwn(w)2 {z^s + z^i) 

x[X+{co,uj)+X_{oj,uj)]Jx^. (44) 

Here, J^^ are the components of the 2-by-2 exchange 
matrix. Substituting Eq. (44) back into Eq. (41), one 



finds the biphoton amplitude for Type II conversion is 
given by 



hE^L 



^i^{2zp+Zds+Zdi) 



Jxn 



A*(cj)^ 



sin A A: f 
A/c# 



[X+(cj,cj)+X_(cj,cj)], 



(45) 
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with the corresponding coincidence count rate 



X I + X_(cL;,a;)] 1^ 



|A(c.)|4 



sin Afcf 



(46) 



IV. DISCUSSION 



Before discussing the details of the count rates in 
Eqs. (39) and (46), there are a few general features worth 
commenting on. Both Type I and Type II rates are 
proportional to the pump intensity, the square of the 
non-linear susceptibility and have a 'sine '-function de- 
pendence on the longitudinal phase matching condition, 
all of which are expected features of PDC rates. Whilst 
the first two aspects are unaffected by absorption, the 
'sinc'-function, in certain cases, can be distorted (see 
Fig.[4|. 




FIG. 4: Dependence of the longitudinal phase matching con- 
dition for vanishing absorption (solid line), non- vanishing ab- 
sorption (small dashed line) and the special case of degen- 
erate, frequency independent absorption (large dashed line). 
/ = I sine [A/c^] e*^'^ 2" The inset shows an expanded view 
of the first minimum. The units on both axes are arbitrary. 

In the presence of absorption, in addition to the ex- 
pected reduction in amplitude one finds a raising of the 
minimum. This is caused by the reduction in intensity 
of one mode compared to the other, which leads to in- 
complete destructive interference. In the special case of 
degenerate conversion and identical absorption at all fre- 
quencies one sees complete destructive interference be- 
cause, in this case, the amplitude of all modes are reduced 
equally. 

The inclusion of absorption adds two new features to 
the count rate. Firstly, there are additional interaction 
terms, such as the nonlinear interaction of the electric 
field with the noise polarization field. Secondly, one ob- 
serves the expected exponential decay of the modes via 



linear absorption. In the following we will consider the ef- 
fects of PDC in beta barium borate, p — BaB20/^ (BBC), 
a common nonlinear crystal that can produce both Type 
I and Type II down conversion depending one the ori- 
entation of the pump beam to the crystal axis. In real- 
ity, BBC is strongly birefringent (c.f. Table [l| see also 
Ref. |23J). However, for our example we neglect this as- 
pect as it does not affect the features that are due to 
absorption; the features we wish to highlight in this arti- 
cle. Thus, we take the refractive index n{ijo) to be equal 
to the values of the ordinary wave, n(co'o), as opposed to 
its value of the extraordinary wave, nioOe)- 

Let us first consider the effect of the extra interac- 
tion terms, which refer to the down- conversion of the 
pump field to one or more noise polarization fields. From 
Eqs. (19) and (20), one observes that the contribution of 



these extra interactions to the rate is contained in the 
A{uj) prefactors as defined in Eq. ( [21] ). For vanishing ab- 
sorption, A{ijo) 1. As shown in Fig. [sj as absorption 
increases, so does A{uj). The increase is due to the extra 



\A{oj) 

5.x 10 









4.x 10-12 


\ \ 








N \ \ 






3.x 10-12 


''^^ \ 




/ y 
















2.x 10-12 



- 1 




2.x 10^ 4.x 10^ 6.x 10^ 8.x 10^ LxlG^^ 



FIG. 5: Relative increase of A(uS) with absorption, n' is the 
imaginary part of the refractive index. The real part of the 
refractive index is taken to be Re\p\ — 1.67. 

conversion pathways available to the pump mode. This 
increase in final-state phase space causes an increase in 
the probability of conversion and hence the rate. This 
kind of transition enhancement effect has already been 
seen in linear systems, for example in the surface-induced 
transition rates in Refs. [191 [24]. For significant absorp- 
tion of around 10% per cm we expect a relative increase 
in the count rate on the order of 5 x 10~^^. 

Superimposed on the small rate increase caused by the 
extra interaction terms is the effect of linear absorption 



A(nm) 


Re\n(ijJo) 


Re\n(uje)\ 


1064 


1.65 


1.54 


532 


1.67 


1.55 


266 


1.75 


1.61 



TABLE I: Table of frequency dependent refractive indices for 
BBO [23,. 



on the propagating modes. Of the two effects, this is 
by far the more dominant. Figure [6] shows the rela- 
tive change in the count rate as absorption is increased. 
Notice that the Type I and Type II rates differ. This dif- 
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FIG. 6: Relative change of the coincident count rate with ab- 
sorption for Type I (solid line) and Type II (dashed line) con- 
version for a crystal length of L = 2mm. Here, Re[n{2uj)] — 
1.75, i?e[n(cj)] = 1.67, cj = 3.54 x 10^^s-\ The imaginary 
part of the refractive index n" is taken to be frequency inde- 
pendent. Rq is the rate for vanishing absorption. 



ference is due to the interference between the pump beam 
and the outgoing modes, and between forward propagat- 
ing modes and those which have been reflected from the 
crystal interfaces [c.f. Eqs. (|33|, (|34|, (|3^ (|46|]. In 
fact, as one changes the crystal length the relative rates 
of Type I and Type II change, with some lengths result- 
ing in a larger Type II rate compared to Type I. One also 
observes the appearance of beat frequencies in the count 
rate as the crystal length is changed (see Fig. [7|. These 
are also caused by interference of the modes within the 
crystal. Note also that a longer crystal length leads to 
a larger distance over which linear absorption can act, 
hence the general reduction in the count rate with in- 
creasing crystal length. 

Let us comment on a few observations. Firstly, owing 
to the dyadic structure of the Green tensor, the detected 
modes and the modes created at the interaction point 
decouple. Thus, the count rate of a specific polarization 
state is a function of both possible polarization states 
[c.f. Eqs. (23) and (40)]. This highlights a fundamen- 
tal principle of quantum mechanics that, until measured, 
the signal and idler modes propagate as a superposition 
of both possible polarization states. Another observation 
is that the contributions to the count rate from each of 



the detected polarization states are equal [c.f. Eqs. (38) 



and (45)]. It is easy to show that, despite absorption, 
this is a maximally entangled state with respect to po- 
larization. In fact, since the outgoing modes propagate 
as both polarization states, even with a more compli- 
cated material, (e.g birefringence, frequency dependent 
absorption, etc.) the detected state will still be max- 
imally entangled. Similar arguments can be used show 
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(a) Type I Conversion 
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(b)Type II Conversion 

FIG. 7: The interference of the modes within the crystal for 
Type I and Type II down conversion. The beat frequencies oc- 
cur with changing crystal length as the reflected mode moves 
in and out of phase with the propagating mode. Parameters 
as above. 



that the Hong-Ou-Mandel interference pattern from two- 
photon interference is unaffected and it is still possible to 
get maximal visibility despite absorption. It is important 
to note that this analysis has described a experimental 
setup that post-selects the two photon state. It is possi- 
ble that during the propagation through the crystal one 
photon of the down-converted pair is absorbed and lost 
completely. However, these states are not detected in 
coincident counting experiments and hence do not affect 
the results outlined above. 



V. SUMMARY 

In this article we have investigated the coincident count 
rates for Type I and Type II parametric down- conversion, 
subject to absorption. Despite initial assumptions to the 
contrary, we found that the increase in rate that is caused 
by the extra nonlinear interactions of the electric field 
with the noise polarization field are negligible compared 
to the effect of linear absorption. Thus we have con- 
clusively shown that nonlinear absorption and nonlinear 
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noise interactions can be safely neglected in all current 
experimental work. We also found that the rates for Type 
I and Type II conversion differ. Aside from the possibil- 
ity of differing nonlinear susceptibilities the two types 
of conversion, they also differ owing to the different in- 
terference processes that occur between the waves that 
scattered from the crystal boundaries. 

We also observe that, in spite of absorption, the en- 
tanglement and two-photon interference properties of the 
photon pairs are unchanged by both the extra nonlinear 
noise interactions and by linear absorption if one assumes 
that both photons are detected. However, absorption will 
play an important role in single photon counting exper- 
iments where one only detects one of the photons in the 
pair. Unlike in coincident counting experiments, where 
the loss of one photon means that the photon pair is un- 
detected, in the case of single photon counting these par- 
tially absorbed photon pairs will contribute significantly 
to the count rate. Hence, in single photon counting ex- 
periments both linear absorption and nonlinear noise in- 
teractions will play a more important role. 
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to be in region 2 and region 3, respectively, one always 
has Zd > za- Hence, in this case the positive sign is valid. 

The 2;-dependent factors are given by 

X [e-'^^'^ ^rl\uj)e'^^^'^^^^]M^{uj). (A3) 
with a TE^ TM. The relevant Fresnel coefficients are 



21/23 



(to) 



kz - Qz 21/23. _kz- e{uj)q 



kz + Qz' kz + e{oj)qz ' 

t¥Ei^) = T:^, tfMi^) = T-^^- (A4) 



w ith the m agnitude of the wave vector in vacuum Qz = 
^Jq^ — k\ with q = uj/c. The function 



Mte/Tm{^) - - ^TE/TAli^yTE/TMi 



accounts for multiple reflections inside the crystal. 
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Appendix B: Commutation relations 



Appendix A: The Green tensor for planar 
multilayered media 

We consider a dielectric medium that is layered in the 
2;-direction and consists of three regions (vacuum-crystal- 
vacuum). Each region is infinitely extended in the (x, y)- 
plane. The interfaces between the crystal and the vac- 
uum are located at z = ±L/2. The scattering Green 
tensor that describes transmission through the output 
face of a crystal at z = +L/2 is given by [25j 



Gc/3(rd,rA,cj) 



2(2^)^ 



/ 



d^k± 1 



kl 



[Ma{k^)Mp{-k^)F^%{za, za) 
- N^{k^)Np{-k^)F^\t{za.ZA% (Al) 



with TE and TM vector wave functions 



M(ki) = i(ki X 5), iV(ki) 



-(kj± X k^ X S) 



(A2) 

Here = (x,^, 0) and = {kx^ky^O) are vectors 
restricted to the (x, ^)-plane, k± = |k^| and kz = 
yjk^ — k\ with k = n{rATOj)uj/c being the magnitude 
of the wave vector in the crystal. The vector kj± is de- 
fined as kj± = k^ ± /c^i where the upper (lower) sign is 
taken for Zd > za {zd < za)- As za and Zd are assumed 



The commutation relations for the macroscopic fields 
can be found from the bosonic expansions of the electric 
and noise polarization fields, Eq. Q and Eq. Q, respec- 
tively. On application of Eq. ^ the commutators for the 
macroscopic fields are found to be 



£:«(r,a;),^^(r^cJO 



^«(r,a;),P^^^(r^cJO 



7T€oC 



2 ImGa/3(r,r',cj)(5c^c^/, 

(Bl) 



TTC^ 



Pjv,a(r,w),£;t(r',w') 



e"{r', oj)Gai3{r, r', oj)S^^' , 
(B2) 

— ^s''(r, a;)G;^(r, r',uj)S^^^ , 

(B3) 

PAr,a(r,cj),Pj^^(r^cj')j = -Eoe" {r,uj)S{r - r')Sa(3^u;oj' - 

(B4) 



Appendix C: Products of vector wave functions 

The contraction of the TE and TM vector wave 
functions with the susceptibility for a Type I process, 
da/si^si^i) = d{(jUs,uJi)[xx + yy]^i^i Icads to the follow- 
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mg expressions 

Ms,a{-'k±)[xx + yy]^^Mi,p{k_i_) = kl, 

M,,„{--k_i_)[xx + yy]^pNi,l3{k_i_) = 0, 



(CI) 



(C2) 



iV«,a(-k±) [xx + yy]^^Ni,^{k^) = kl 



(C4) 



Similarly, for the susceptibility for a Type II process, 
da/3{(^s,(^i) = d{ujs,uji)[xy -\-yx]^^, one finds 



Ms^a{-^±)[xy -\- yx] Mi^^{k^) = -2k^ky, 



(C5) 



M,,,(-kx) [xy + yx] ,^iV,,^(kx) = i{kl - k^ (^^) 



f k k 

^^s,a(-k_L)[x^ + p]^^iVi,^(k_L) = 2k^ky ( 



(C6) 

kg 
(C7) 



(C8) 



The (dyadic) tensor products of the uncontracted vector 
wave functions are 



M^,„(k_L)Mi,^(-k_L) 



/ kl 



kxky 



kx ky k^ 

V 0, 



(C9) 



(C3) M,,«(ki)A^i,/3(-ki) 

2 / '^kxkykz^i 



ki 



zkykj^ ^ 



ikykz^i 



/ 



(CIO) 



ikxkykz^s '^k'^kz^i 
- I ikykz^s ~'^kxkykz^i | , 
iky k I 2/kx k I 



(Cll) 



A^.,«(k±)Ari,;3(-ki) = 

i^x z,s'^z,i '^x'^y'^z^s'^z^i i^x i^z ,s _\_ 

■ ■ 

rvxi^yi^z^si^z^i i^yi^z^si^z^i '^y'^z,s'^_\_ 



ka ki 



\ —h h —h h 

\ '^x'^z,i'^_\_ '^y'^z,i'^j_ '^_\_ 



(C12) 
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